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Abstract 

We present the temperature-dependent elastic constants of two ductile rare-earth intermetallic compounds 
YAg and YCu with CsCl-type B2 structure by using a first-principles approach. The elastic moduli as a 
function of temperature are predicted from the combination of static volume-dependent elastic constants 
obtained by the first-principles total-energy method with density functional theory and the thermal expan- 
sion obtained by the first-principles phonon calculations with density-functional perturbation theory. The 
comparison between our calculated results and the available experimental data for Ag and Cu provides 
good agreements. In the calculated temperature - lOOO/T, the elastic constants of YAg and YCu follow a 
normal behavior with temperature that those decrease with increasing temperature, and satisfy the stability 
conditions for B2 structures. The Cauchy pressure for YAg and YCu as a function of temperature is also 
discussed, and our results mean that YAg and YCu become more ductile while increasing temperature. 
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1 Introduction 

The rare-earth intermetallic compounds typically possess high strength and stiffness, low specific weight corrosion re- 
sistance, and hot strength superior to ordinary metals |[T]|2l. However, their room temperature brittleness and poor fracture 
toughness severely restrict their applications. Recently, a large family of rare-earth intermetallics has been discovered to 
be ductile and tough at room temperature |3 1, in some cases exceeding 20% ductility in tension. These alloys have the B2 
structure in space group Pm3m (number 221) and chemical formula RM, where R denotes a rare-earth element, and M de- 
notes a late transition metal or an early p-element. Many theoretical and experimental works P1l5ll6ll7l[8ll9l [T0lfTTIfT2l[T3]| 
have been done to understand the mechanism of high ductility at room temperature, such as elasticity, defect properties, 
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phase stability, dislocations, stacking fault, electronic structure and density of states (DOS), etc. However, the properties 
and mechanism have not been understood completely hitherto. 

To get a better understanding of the anomalous ductility of RM intermetallics, more fundamental investigations of 
temperature-dependent properties such as thermal expansion, heat capacity, and thermoelasticity, etc., are obviously re- 
quired. In general, elastic properties of a solid are very important because they are closely associated with various 
fundamental soUd-state properties such as interatomic potentials, equation of state, and phonon spectra. The temperature 
dependence of the elastic constants of a material is important for predicting and understanding the mechanical strength, 
stability, and phase transitions of a material lfT4l . However, there are few literatures focused on the temperature-dependent 
moduli of RM intermetallics. The first-principles methods based on density functional theory (DFT) can compute the 
single-crystal elastic moduli at zero-temperature accurately, but treating the corresponding temperature dependence of 
elastic constants is still a formidable challenge ifTSl . More recently, Wang et. al. lfT6l[r7l proposed a simple quasistatic 
approach to determine the elastic constants at finite temperatures by using first-principles calculations, and the excellent 
agreement between the calculated results and experimental data is found. Their approach is based on the fact that the tem- 
perature dependence of elastic constants mainly results from volume expansion with increasing temperature lfT4l [TSl [T9l . 

In the present work, we apply the first-principles quasistatic approach to investigate the temperature-dependent elas- 
tic constants for the novel intermetallics YAg and YCu. The static volume-dependent elastic constants are obtained 
by using the DFT total-energy calculations combined with the method of homogeneous deformation. A first-principles 
quasiharmonic approach Il20l is applied to predict the phonon free energy and volume versus temperature by using density- 
functional perturbation theory (DFPT). To benchmark the reliability results of the presented method, we have also calcu- 
lated the temperature-dependent elastic constants for Ag and Cu, and the comparison between our predicted results and 
the available experimental data provides good agreements. In addition, the ductility of YAg and YCu as a function of 
temperature has been investigated by using the temperature-dependent elastic modulus. 

2 Theory 

In this section, we present the necessary theory for the calculation of the temperature dependence of cubic crystal 
elastic moduli in B2 structures. For programming the calculation of the elastic constants, we consider a procedure of 
homogeneous deformation of a bulk-crystal. With the displacement given as m, = x, - a, (/ -1,2, 3), between the initial 
configuration a, and the strained configuration jc,, the deformation applied to the crystal is described as the deformation 
gradient matrix Jjj - dxjdaj. Then, we may define the Lagrangian strain T]ij - \(JikJkj - Sij) |2 1 1. Thus, the isothermal 
elastic constants Cjjj^i are defined as the derivative with respect to the Lagrangian strain on the Helmholtz free energy 
Eq.([8]) at constant temperature ll22l . 

-1^' (1) 
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where V - V{T) is the equilibrium volume at temperature T, and 77' indicates that all other stains are held fixed. This 
is the general way to compute the finite temperature elastic constants in principle. However, this procedure is quite 
cumbersome at the present since it involves calculation of the second derivatives of F^w,. For cubic materials, there are 
there independent elastic constants C^-^, Cjj, and Cj^ (in Voigt notation) which describe the elastic behavior completely. 
A more convenient set for computations are Cj^ and two linear combinations and jj.^ . The bulk modulus 



= (C[i + 2Cf2)/3, 

is the resistance to deformation by a uniform hydrostatic pressure; the shear modulus 

fi'' = (C[i - C[2)/2, 



(2) 



(3) 



is the resistance to shear deformation across the (110) plane in the (110) direction, and C44 is the resistance to shear 
deformation across the (100) plane in the (010) direction. The bulk modulus B is determined from the Vinet equation of 
state ||23l . The shear modulus fi is calculated from volume-conserving orthorhombic strain 



(60 
6 
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and the Helmholtz free energy related to this strain is 



F(V, 6) = F(V, 0) + 6f/V6^ + 0(6^), 



(4) 



(5) 



where F(V, 0) is the free energy of the unstrained structure. We use a volume-conserving tetragonal strain to determine 
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which leads to the free energy change 



F{V,6) = F(y,0) + 2Cl4V6^ + 0(6'^), 



(6) 



(7) 



In the usual way and following Moriarty et al L24il . the Helmholtz free energy F for a crystal at temperature T and 
volume V can be separated as 

F(V, T) = £„(y) + FM T) + Fei(y, T), (8) 

where Eq{V) is the zero-temperature total energy of the electronic ground state, Fvib(V, T) is the vibrational free energy 
which comes from the phonon contribution, and where Fgi represents the thermal electronic contribution to free energy 
from finite temperature. Here V is taken as the volume per unit cell and F is the free energy per unit cell. The total 
hydrostatic pressure in the crystal can be calculated by solving 

dF{V, T) 



P{V,T) = -- 



dV 



(9) 



with corresponding relations for its three components. Other thermodynamic functions can be obtained from F{V, T) 
similarly, e.g., entropy S - -(dF/dT)v, internal E = F + TS , and the Gibbs free energy G - F + PV. 

In usual way, the vibrational free energy Fyib{V, T) includes both quasiharmonic and anharmonic components. The an- 
harmonic effects which represent the phonons interaction and electron-phonon coupling are so weak and can be neglected 
ifTS l. In this work, we only consider the remaining quasiharmonic phonon free energy which can be written as 



(10) 



qA L 

where represents the Boltzmann constant, H is the reduced Planck constant, and ojq^ represents the frequency of the 
Ath phonon branch at wave vector q. 

With finite temperature DFT calculations 1251 . the thermal electronic contribution to free energy F^i is obtained from 
the energy and entropy contribution and given by 

Fel = £el - r5el, (11) 

with the bare electi^onic entropy Sei 

SdV,T)^-ks n(s,V)[f(s)\nf(s) + (l-f(s))\n(l-f(e))]de, (12) 
Jo 



and the thermal electronic energy Eei 



E,i(V,T) 



I n{E,V)f{e)edE - \ n(s,V)£de, (13) 
Jo Jo 



where n(s, V), /(e), and Sf represent the electronic density of state (DOS), the Fermi-Dirac distribution, and the Fermi 
energy, respectively. 

In present work, we employ a first-principles quasistatic approach, which was developed by Wang et al. lfT6l flTl . 
to calculate the temperature-dependent elastic constants Cij(T). Following the quasistatic approximation, the change of 
elastic properties at elevated temperatures is mainly caused by volume change due to thermal expansion, and the contribu- 
tions of Fvib and Fe\ to these second derivations can be neglected. In fact, the contributions due to the kinetic energy and 
the fluctuation of microscopic stress tensor are small and can be ignored reasonably |26l. The thermal expansion and the 
equilibrium volume V at T are obtained by using the quasiharmonic approximation in our calculation. A similar procedure 
was employed by Kadas et al EtI . though they calculated thermal expansion by using Debye-type model. We obtain the 
temperature dependence of isothermal elastic constants Cjj(T) by the application of the following three-step procedure. 
The first step is calculating the thermal expansion and the equihbrium volume V(T) at T by using the first-principles 
quasiharmonic approach. In this procedure, one can compute volume-dependence of phonon density of states at OK in a 
set of volume points, and the predicted equilibrium volume V{T) at T ^ OK (or inversely T(V) relation) is determined 
from fitting to the Vinet equation of state ll23l by minimizing free energy F with respect to V, while isothermal bulk 



modulus as a function of temperature T is also obtained. In the second step, we obtain the volume-dependent elastic 
constants Cfj{V) atT - OK as the second derivatives of the Helmholtz free energies with respect to strain tensor by using 
the energy-strain relation based on Eqs.(2]|7i. In the third step, the calculated elastic constants from the second step at 



the volume V{T) are approximated as those at finite temperatures, i.e., Cjj{T) = Cjj(T(V)). There is ample experimental 
evidence 1281 [T9l |29]| to support this approximation in which the temperature dependence of elastic constants are solely 
caused by thermal expansion. 

In order to compare with experimental data, the isothermal elastic moduli Cfj must also be transformed to the adiabatic 
elastic moduU C?. by the following relation Il30l : 



-'7 U Cv ' ^ 
where Cv is the specific heat at constant volume and 



. T TV 

C i - Cl + —A,Aj (14) 



Ai^Y.'^'Cl (15) 



with the Unear thermal expansion tensor ■ For cubic crystals, Eq. ( 14 1 simplifies to 



Cl = C^, (16) 

Cfi=C[i+A, (17) 

Cf2 = C[2 + A, (18) 

where ifTSl 

TV t2 

A = —a^B^\ (19) 

with the volume thermal expansion coefficient a - V^^dV IT\p. It is worth noted that the correction A increases with 
temperature. 

3 Computational details 

In present work, the static energy and the thermal electronic contribution to the Helmhotz free energy were computed 
by using the first-principles calculations in the framework of the density-functional theory (DFT). We employed the plane- 
wave basis projector augmented wave (PAW) method ll3n[32l within the generalized gradient approximation (GGA) in the 
Perdew-Burke-Ernzerhof (PBE) Il33ll34l exchange-correlation functional as implemented in the VASP code 1351 l36l [37l . 
The radial cutoffs of the PAW potentials of Ag, Cu, and Y were 1.50, 1.31 and 1.81 A, respectively. The Ad and 5s 
electrons for Ag, the 3d and As electrons for Cu, and the As, Ap, Ad and 5i electrons for Y were treated as valence and 
the remaining electrons were kept frozen. The Brillouin zones (BZ) of the unit cells are represented by Monkhorst-Pack 
special A: -point scheme ll38l . Since high accuracy is needed to evaluate the elastic constants, the convergence of strain 



energies with respect to the Brillouin zone integration was carefully checked by repeating the calculations for 21x21x21 
and 25 X 25 X 25 grid meshes at equilibrium volume at OK, and we found at most 0.5GPa difference both for C44 and fi. 
Hence, we used 21 x21 x21 in the full Brillouin zone giving 726 irreducible k-points. In addition, we used a high plane- 
wave energy cutoff of 600eV which is sufficient to calculate the elastic moduli accurately. The thermal electronic energies 
and entropies are evaluated by using one-dimensional integrations from the self-consistent DFT calculations of electronic 



DOS using FD smearing as shown in Eqs. ( 12 1 and ( 13 1. Calculations of the static strain-energy were performed by the 
tetrahedron method with Bloc hi corrections 1391 , which gives a good account for total energy. 

The vibrational free energy were obtained from the first-principles phonon calculations by using PHONOPY Il20l 
l40l l4n which can support VASP interface to calculate force constants directly in the framework of density-functional 
perturbation theory (DFPT) B2ll . Phonon calculations were performed by the supercell approach. Since the chosen 
supercell size strongly influences on the thermal properties, we compare the vibrational free energies of 3 x 3 x 3 supercell 
with those of 5 x 5 x 5 supercell at 300K and lOOOK, and find that the energy fluctuations between 3x3x3 and 5x5x5 
supercells are less than 0.01%. Hence, we chose the 3x3x3 supercell with 54 atoms to calculate phonon dispersions. 
We carried out DFPT calculations on this 54 atoms supercell using PBE-GGA exchange-correlations effects and 7x7x7 
k-point grid meshes for BZ integrations. In order to deal with the possible convergence problems for metals, a smearing 
technique is employed by using the Methfessel-Paxton scheme Il43ll . with a smearing with of 0.05eV. 

We have calculated Helmholtz free energy [the right-hand side of Eq.(|8]l] at temperature points with a step of IK from 
to lOOOK at 15 volume points. At each temperature point, the equilibrium volume V and isothermal bulk moduli are 
obtained by minimizing free-energy with respect to V from fitting the integral form of the Vinet equation of state (EOS) 
Il23l . Then, the volume thermal expansion coefficient was obtained by numerical differentiation for dV/dT. 

4 Results and discussion 

Table jTjgives our calculated lattice constants and isentropic elastic constants Cfj atT - OK and T - 300K at ambient 
pressure {P - OGPa) in comparison with the results from the previous calculations [] (T = OK) and experimental data at 
ambient conditions ||4l(f = OGPa, T = 300K). In all cases, the comparison is quite agreeable. At T = OK our calculated 
results show a very good agreement with the other theoretical data. For YCu, the present results at 300K are very closed 
to the experimental values obtained by Morris et al [4| at ambient conditions. The largest discrepancy between theory and 
experiment is the value for C44, where the difference is approximately 10% for YAg. 

For the temperature dependence of the isentropic elastic moduli C? , we present in Figures Tjj4 1 



our calculated in the 



range 0-lOOOK at ambient pressure. Figures [T] and |2] show our findings for benchmark metals Ag and Cu respectively, 
accompanied by available experimental data taken from ultrasonic measurements ll44l l45l l46l . Figures |3] and |4] give our 
prediction for the unknown values of temperature-dependent C^j for YAg and YCu, respectively. Through the calculated 
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values of isentropic elastic constants for Ag and Cu in comparison with the values measured by experiments, we get 
a useful test of the accuracy of the method and the precision of our calculations of temperature dependence of elastic 
constants for YAg and YCu. The same DFT method had also been employed to calculate the elastic constants of Cu as a 
function of temperature by Wang et al 1 16|, and our results show no discrepancy from their values. The overall observation 
is that all the calculated values of C^p C'^j' ^i^^ decrease with increasing temperature, since thermal expansion may 
soften the elastic moduli at high T. At higher temperature, we also find that trend of Cfj approach linearity. Considering our 
calculated results both for YAg and YCu, we find the values of C^j decrease to the largest extant in the whole temperature 
range 0-lOOOK, and those of decrease least. Cfj, Cfj, and C^^ for YAg decrease by 12.8GPa, 7.3GPa, and 9.6GPa, 
and those for YCu decrease by 13.9GPa, 5.9GPa, and 10.3GPa, respectively. The requirement of mechanical stability in 
a cubic crystal leads to the following restrictions on the elastic constants, Cu - Cn > 0, Cu > 0, and C44 > BTI . Our 
results of the elastic constants as shown in Figure [3] and |4] satisfy these stability conditions in the range 0-lOOOK. 

Next we discuss the temperature dependence of Cauchy pressure which was suggested that it could be used to describe 
the angular character of atomic bonding in metals and compounds by Pettifor B8l . The isentropic Cauchy pressure is 
defined as 

P'c = Cn - Cl- (20) 

If the value of Cauchy pressure is more positive and bigger, the bonding of material is more metallic in character. In 
contrast, negative Cauchy pressure requires an directional character and low mobility in the bonding. Generally, for 
ductile materials such as metals Ag and Cu, the Cauchy pressures have positive values, while for brittle semiconductors 
such as Si, the Cauchy pressure is negative. The YAg and YCu are a new class of highly ordered and ductile intermetallics, 
especially for YAg in some cases exceeding 20% ductility in tension |4|. In table[T] the present results of Cauchy pressure 
T = OK and T = 300 K have positive values for YAg and YCu, and agree well with the previous theoretical results at 
T = OK and the experimental results at room temperature (about 300K). Figure |5]illustrates that the Cauchy pressure for 
YAg and YCu as a function of temperature. In the range of temperature 0-lOOOK, the Cauchy pressure of YAg is always 
greater than that of YCu, and those both for YAg and YCu increase with elevating temperature. Our calculated results 
demonstrate that high temperature may soften the directional character, and YAg and YCu will become more ductile while 
the temperature is increased. 

5 Conclusions 

In this work, we present the temperature-dependent elastic constants of two ductile rare-earth intermetallic compounds 
YAg and YCu with CsCl-type B2 structure by using a first-principles quasistatic approach, in which the static elastic 
constants as function of volume are determined by the first-principles DFT total-energy calculations within the framework 
of the method of homogeneous deformation and the volume-temperature relation determined by first-principles phonon 
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calculations based on DFPT with quasiharmonic approach. To benchmark the reliabihty results of the presented method, 
the comparison between our calculated results for Ag and Cu with available experimental data has been performed and 
shows good agreements. In the whole temperature range 0-lOOOK, the elastic constants follow a normal behavior with 
temperature that those decrease with increasing temperature, and satisfy the stability conditions for B2 structures. The 
values of C^j decrease to the largest extant, and those of decrease least. The temperature dependence of Cauchy 
pressure is also discussed. The Cauchy pressures both for YAg and YCu increase with elevating temperature, and that of 
YAg is always higher than that of YCu. Our results mean that increasing temperature may improve ductiUty of YAg and 
YCu. 
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Table 1: The present calculated lattice constants a (A), isentropic elastic constants C^. (GPa), and Cauchy 
pressure P^(GPa) for YAg and YCu at T = 0^ and T = 300K compared to previous computed results and 
experimental data. Note that all the previous calculated results are T = OK values, while the experiments are 
room temperature values about T = 300K. 

Present calculation Previous calculations Experiment 



a 

'-44 
pS 

YCu 



3.645", 3.656" 
98.9", 95.7" 
52.4^51.4" 
35.5", 32.7^ 
16.9", 18.7^ 



3.634^, 
105.0^ 
50.0^ 
37.0^ 
13.0^ 



3.627'' 

102.5^ 

56.5'' 

37.8'' 

18.7'' 



3.619^ 
102.4'^ 
54.0^ 

37.2'^ 
16.8^ 



a 

pS 

_Lc_ 



3.484", 3.492" 
117.4", 114.6'' 
45.9", 45.6" 
36.3", 33.4" 
9.6", 12.2" 



3.477'', 
116.0'^ 
47.0^, 
35.0^ 
12.0'' 



3.472" 

113.6'' 

48.4'' 

36.8'' 

11.6^^ 



3.477" 

113.4*^ 

48.4'^ 

32.3*^ 

15.9" 



" This work at T ^ OK; 
" This work at 7 ^ 300i<:; 

" Ref.[4J, obtained from first-principles calculations at T 
Ref.|7|, obtained from first-principles calculations at T 
" Ref.|i4J, obtained from experiments at T = 300 A'. 
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Figure 1: For single-crystal Ag, the present calculated isentropic elastic constants and corresponding exper- 
imental data. The solid, dashed, and dashed-dotted curves denote the present values of C^^, and C^^, 
respectively. The open symbols denote the corresponding values of ultrasonic measurements by Neighbours 
and Alers [[44 J . 



10 



CO :=. 

o 



140 
12 
100 



60 
40 
20 





200 



400 600 
Temperature [T] 



800 



1000 



Figure 2: For single-crystal Cu, the present calculated isentropic elastic constants and corresponding exper- 
imental data. The solid, dashed, and dashed-dotted curves denote the present values of C^^, C^^ ^44' 
respectively. The open symbols denote the corresponding values of ultrasonic measurements by Overton et al 
Il45l . and the filled symbols display those from the measurements of Chang and Himmel B6l . 




Figure 3: The predicted isentropic elastic constants as a function of temperature for YAg. The solid, dashed, 
and dashed-dotted curves denote the present values of C^^, C^^ ^^'^ C^^, respectively. Note that Cu > Cu, 
Cii > 0, and C44 > at all temperatures, which is consistent with stability of B2 lattice. 
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Figure 4: The predicted isentropic elastic constants as a function of temperature for YCu. The solid, dashed, 
and dashed-dotted curves denote the present values of C'^j, C^^ ^^'^ respectively. Note that Cn > C12, 
Cii > 0, and C44 > at all temperatures, which is consistent with stability of B2 lattice. 
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Figure 5: The predicted isentropic Cauchy Pressure defined by = - as a function of temperature. 
The solid and dashed curves denote the values for YAg and YCu, respectively. 
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